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Abstract 

We classify dark energy models in a plane of observables that correspond to the common 
parameterization of a non-constant equation of state, w(a) = wo + w a (l — a), where a is the 
scale factor of the universe. The models fall into four classes and only two of these classes have 
a region of overlap in the observable plane. We perform a joint analysis of all Type la supernova 
(SNIa) data compiled by the High-Z SN Search Team (HZT) and the Supernova Legacy Survey 
(SNLS) and find that no class of models is excluded by current SNIa data. However, an analysis 
of large scale structure, Lya forest and bias constraints from SDSS, the Gold SNIa data and 
WMAP data indicates that non-phantom barotropic models with a positive sound speed are 
excluded at the 95% C. L. 
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Abundant cosmological data indicate that the expansion of the universe has changed from a 
decelerating phase to an accelerating phase in the last few billion years. This is generally attributed 
to the recent dominance of an energy component with negative pressure called dark energy 1 . For 
reviews see Ref. ■2]. 

Among candidates for dark energy are the cosmological constant A proposed by Einstein and 
a dynamical scalar field such as quintessence 0110. The cosmological constant arises in particle 
physics as vacuum energy with constant energy density p, constant pressure p and equation of 
state w = p/p = —1. A cosmological constant is troublesome from the particle physics standpoint 
because the cosmologically measured value of p 1 / 4 is found to be about ( 0.7p c ) 1/l4 — 0.002 eV, where 
p c is the critical density for which the universe is flat. Quantum field theory, however, suggests a 
value at least 15 orders of magnitude larger under the assumption that the Standard Model is an 
effective theory valid below 1 TeV. For the vacuum energy to be renormalized to 0.002 eV, a fine 
tuning in the energy density at the level 60 decimal places is required. See Ref. 0 for a review of 
the cosmological constant problem. Compounding this, is the cosmic coincidence problem which 
seeks an explanation for why the dark matter and vacuum energy densities are comparable today 
although their ratio scales as 1/a 3 , where a = 1/(1 + z) in terms of redshift z |2|. Quintessence 
is a not a solution, but provides a phenomenological explanation for these problems. It is a time- 
varying held with w(a ) > — 1 and is usually represented by a very light scalar held rolling down 
a potential. The potential and held values are chosen so that the energy density in quintessence 
dominates the dark matter density only recently and takes on the measured value. This picture 
has been generalized to schemes with w < — 1 called phantom dark energy and to dark energy as a 
barotropic fluid with w = f(p)/p. In what follows, we consider phantom models, barotropic fluid 
models of dark energy and two classes of quintessence models. 

Recently, in a sequence of papers 0 0 0, models of dark energy have been divided into 
categories depending on their equation of state, w, and its derivative with respect to the logarithm 
of the scale factor, w' = dfe since H 3 = din(a) ’ ^ time derivative of w in units of Hubble 

time. It was pointed out that different classes of models evolve in different regions of the w — w' 
plane 2 . 

Building on this work, we classify models of dark energy in a parameter space of observables. 
It has been argued in Ref. EH, that given the sensitivities of future cosmological experiments, 
one can realistically expect to constrain at most two parameters related to the equation of state. 

1 Modifications to the Friedman-Robertson-Walker equation can also provide an explanation for the current accel¬ 
erated expansion such as in Cardassian models [TJ which do not have exotic forms of energy or a vacuum contribution. 
We do not consider such models here. 

2 Equivalently, different dark energy models evolve in different regions of statefinder planes ESI; the statefinder 

diagnostic, defined in terms of the second and third time derivatives of a, is directly related to w and w'. 
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Specifically, we use the standard parameterization of dark energy evolution in terms of the equation 
of state and its derivative with respect to z at the present time m 


% 

w(z) = Wo + W a zr— (1) 

-L \ Z 

with 

dw , 

Wa = lte' z=0 = ~ 2W ^ =1 ' ^ 

Note the relationship between w a and the value of w' at z = 1, the epoch when the dark energy 
contribution to the expansion rate is expected to increase in importance 3 . We classify models in 
the wo — w a plane. The primary advantage of this parameterization is that since the dark energy 
density evolves according to 

p(z) = /9(0)(1 + z)3( l+w 0 +w a ) e -3waj^ ) ( 3 ) 

it is well-behaved from high redshifts until today. Note that for dark energy to be subdominant at 
early times, wq + w a < 0. We do not impose the latter restriction since the parameterization of 
Eq. Q, though well-behaved at high redshift, may be inapplicable there. 

The previously popular linear parameterization w(z) = wq + w\z m diverges at high redshift 
and consequently yields artificially strong constraints on w± in analyses involving data at high 
redshifts such as CMB data. Since Eq. 0 is often used in the analysis of data, our classification 
will enable a direct comparison between the results of such analyses and the theoretical space of 
models. 

The categorization we develop is analogous to that of inflationary models where models popu¬ 
late regions in a plane defined by the spectral index of scalar perturbations, n s , and the relative 
normalization of the tensor and scalar spectra, R or in a plane of horizon-flow parameters Em 

Dark energy models: 

Here, we briefly review the models under consideration, define the regions they occupy in the 
w — w r plane and classify them in the wo — w a plane. 

We first consider models which obey the null energy condition, w > —1, and then consider 
phantom models with w < — 1. 

Thawing models 

The nomenclature, “thawing” models, was coined in Ref. [Jj to describe a scalar field whose 
equation of state increases (thaws out) from w ~ —1 as the scalar rolls down towards the minimum 

3 In Refs. lUH 1141 . the transition from matter domination to dark energy domination is found to occur at z ~ 0.5. 
For a recent discussion of the low redshift evolution of the dark energy density see Ref. EH- 
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Figure 1: Classification of the four types of models in the (wQ,w a ) plane. All regions depicted by 
a light yellow shaded region correspond to phantom models with w < — 1 today or recently. 


of its potential. Potentials of the form cj) n (with n > 0) and e ^ are typical of these models. Particle- 
physics models involving axions ESI, pseudo Nambu-Goldstone bosons j|4j|, moduli or dilatons m 
often have such potentials. They are found to satisfy 7] 

1 + w < w' < 3(1 + w ). (4) 


It is natural to impose this constraint at the epoch when dark energy starts becoming important, 
i.e., z = 1. With rc(l) = wq + w a / 2 and w'\ z= \ = — w a /2 , we find 


- -(1 + Wq) <w a < -(1 + w 0 ). 


(5) 


Thawing models occupy the dark-shaded triangular wedge in Fig. |3 


Cooling models 

In these models, initially, w > — 1 and w decreases as the scalar rolls down the potential, 
which is typically of the form (j)~ n or cj)~ n e^ 2 (with n > 0). The tracker models of Ref. jJT] adopt 
such potentials. Such forms of the potential also arise in models of dynamical supersymmtery 
breaking [22] and supergravity [23]. These models lie in a region of the w — w' plane defined 

by EL IE 9j 

— 3(1 — rc)(l + w)<w'< 0.2rc(l + w). (6) 

By requiring that rc(l) obey the above inequality, we obtain the region with red slant hatches (with 
negative slope) in Fig. ^ As shown in Ref. ^-essence models m with a non-linear kinetic term 
also fall within the class of cooling models. 
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Note that the “freezing” models of Ref. 7j are a special case corresponding to the situation in 
which the potential has a minimum at <f> = oo. For freezing models, 


3ic(l +w)<w'< 0.2u;(l + w). 


(7) 


Barotropic fluids 

A barotropic fluid is one for which the density depends only on the pressure. Supposing that 
dark energy is a barotropic fluid with p = f(p) leads to a class of models that contains as special 
cases the generalized Chaplygin gas HU for which f(p ) = —A/p a with a > — 1 and the orginal 
Chaplygin gas model with a = i EH. While the models of Refs. Earn] were proposed to unify 
dark matter and dark energy, subsequent work addressed the possibility of Chaplygin gas models 
as being models of dark energy alone E3 Barotropic fluid models arise in string theories ESI hi 
which e.g. the Chaplygin gas corresponds to a gas of d-branes in a d + 2 spacetime m- 
Non-phantom barotropic models were found to satisfy jS] 

w' < 3u>(l + w) (8) 

under the assumption that c 2 s = dp/dp > 0 so that perturbation growth is stable. 

The region with green slant hatches (with positive slope) in Fig. depicts barotropic fluids. 
This region partially overlaps with that for cooling models. 

Comparing Eq. m with Eq. dHJ) , we see that freezing models and barotropic fluids have nonover¬ 
lapping regions with a common boundary in the w — w' plane. In the observable wq — w a plane, 
freezing models occupy the region with red slant hatches (with negative slope) below the cross 
hatched region (non-overlapping with the region occupied by barotropic fluid models). See Fig. 

Phantom models 

Finally, we consider phantom models for which w < — 1 today or recently. It was shown in 
Ref. U3 that for the prototypical phantom (he. a ghost with a negative kinetic term) model TIT). . 

3(1 — u;)(l + w) < w' < 3tc(l — u;)(l + w). (9) 

However, we do not impose this inequality since it is very marginally more constraining than 
w < — 1, and since the upper bound on w' was obtained specifically for tracker phantom models. 
Requiring wq < -1 or «;(1) < — 1 we find that phantom models populate the light yellow shaded 
region in Fig. ^ 

The subregion defined by vertical (horizontal) lines corresponds to wo < — 1 and u>(l) > — 1 
(n?o > —1 and ie(l) < —1) he., models for which the equation of state crossed the phantom divide 
line w = — 1 from a higher value to a lower value (lower value to a higher value). The light yellow 
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shaded region without lines corresponds to “pure phantom” models for which the phantom divide 
line has not been crossed recently i.e., wq < —1 and w( 1) < —1. 

Note that not all models that violate the null energy condition require the scalar to be a 
phantom. Demonstrations of the violation of the null energy condition have been made in models 
of vacuum metamorphosis EH, climbing scalar fields EH, and the braneworld EH without the 
introduction of negative energies or negative norm states. 

Interestingly, the quintom model of Ref. EH lies in the region with vertical lines. 

We caution against interpreting this classification as being comprehensive. For example, dark 
energy with an oscillating equation of state EH can rarely be realized in a potential formulation. 
(An example of oscillating dark energy from a quite complicated potential can be found in Ref. [SB]). 
Thus, oscillating dark energy as a class of models falls outside the realm of our classification unless 
the equation of state crosses the phantom divide line as in quintom models EH- Similarly, we have 
not attempted to include neutrino dark energy EH or dark energy with generalized equations of 
state m in the classification. 

When viewed in conjunction with data, our classification is intended to provide guidance as 
to what kinds of schemes are preferred by data. It is necessary that a cosmological constant 
( wo,w a ) = (—1,0), be excluded by data to distinguish between thawing, cooling and phantom 
models because these models have a cosmological constant as a limiting case. This is in contrast 
to barotropic fluids with a positive sound speed. 


Current status: 

Having defined the regions that the four different types of models occupy in the wo — w a plane, 
we now examine if current data are able to discriminate between these classes. 

We perform a joint analysis of SNIa data compiled by the High-Z SN Search Team (HZT) and 
the Supernova Legacy Survey (SNLS). We include the Gold and Silver datasets of Ref. [13, the 
4 SN of Ref. [HI, and the SNLS dataset [S3 assuming a flat universe with 0.2 < < 0.4 and 

0.5 < h < 0.9. We have removed the SN common to both surveys from the SNLS dataset which 
results in a full dataset of 266 SN. Our joint analysis is implemented as follows. The statistical 
significance of a cosmology is determined in terms of x 2 = Xhzt + x|nls, w ihh 

188 (fi° bs - 5log 10 d L (z i -,Q rn ,h,w 0 ,w a ) - Mi) 2 


Xhzt 


xInls 


E 

2=1 

78 

E 

2=1 


o 7 


(nf s — 5 log 10 d L (zi-Q m ,h,w 0 ,w a ) - M 2 Y 
°li + a int 


( 10 ) 

( 11 ) 


where /r ° bs is the distance modulus at redshift Zi, dL is the luminosity distance (in units of 10 pc), 
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Figure 2: 68% C. L. and 95% C. L. allowed regions from a combined analysis of the HZT and SNLS 
SNIa data. The crosshair marks the best-fit point (wo,w a ) = (—1-2,1.5). No class of dark energy 
model is excluded. 

ai is the total uncertainty in the distance modulus, a‘j H is the measurement variance and is 
the intrinsic dispersion of SN absolute magnitudes. M\ and M 2 are nuisance parameters that are 
marginalized over in the fit. They correspond to analysis-dependent global unknown constants in 
the definition of distances. M\ and M 2 are analogous to absolute magnitudes if observational data 
are provided as apparent magnitudes. 

Note that SNLS provides the apparent magnitudes, the stretch factor used to calibrate them, 
and the rest frame color factor that measures host galaxy dust extinction. We have chosen to use 
the derived fi values instead for two reasons: (1) it is not possible to perform an analysis at the 
same level of sophistication for the HZT data since the High-Z SN Search Team provides only the 
distance modulus (i.e., after corrections have been made to the apparent magnitude m(z)), and 
(2) it would be overkill considering that current SN data are not very constraining in the (wo,w a ) 
plane. 

Reference m finds c Ti n t = 0.13 ± 0.02. To be conservative in our combined analysis we choose 
&int = 0.15. We emphasize that our approximate analysis of SNLS data alone yields good agreement 
with that of the more accurate analysis of Ref. m implemented in Ref. m- Gravitational lensing 
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bias could be reduced by use of the flux-averaging procedure of Ref. ga¬ 
in Fig. [2] we display the 68% and 95% C. L. allowed regions obtained in our analysis. We 
immediately conclude that all classes of dark energy models are comfortably allowed by current SN 
data. In Fig. [3 we overlay the results of an analysis of large scale structure, Lya forest and bias 
constraints from SDSS, the Gold SNIa data and WMAP data on the plane of classification. The 
68% and 95% C. L. allowed regions have been adapted from Fig. 9 of Ref. m after accounting for 
the fact that our w a corresponds to —w\ in Ref. :32j. 

It appears that current data exclude non-phantom barotropic models with a positive sound speed 
at the 95% C. L. Also, among cooling models, freezing models appear to be favored. However, it 
may be too soon to draw these conclusions with confidence because the analysis of several correlated 
datasets has led to these results. It will be more convincing if the same conclusions can be made 
from an analysis of fewer distinct but larger datasets. For a recent assessment of the abilities of 
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Figure 3: 68% C. L. and 95% C. L. allowed regions from an analysis of large 
forest and bias constraints from SDSS, the Gold SNIa data and WMAP data. 
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